Modulatory effect of Gracilaria gracilis on European seabass gut microbiota community and its functionality

Seaweeds are an important source of nutrients and bioactive compounds and have a high potential as health boosters in aquaculture. This study evaluated the effect of dietary inclusion of Gracilaria gracilis biomass or its extract on the European seabass (Dicentrarchus labrax) gut microbial community. Juvenile fish were fed a commercial-like diet with 2.5% or 5% seaweed biomass or 0.35% seaweed extract for 47 days. The gut microbiome was assessed by 16S rRNA amplicon sequencing, and its diversity was not altered by the seaweed supplementation. However, a reduction in Proteobacteria abundance was observed. Random forest analysis highlighted the genera Photobacterium, Staphylococcus, Acinetobacter, Micrococcus and Sphingomonas, and their abundances were reduced when fish were fed diets with algae. SparCC correlation network analysis suggested several mutualistic and other antagonistic relationships that could be related to the predicted altered functions. These pathways were mainly related to the metabolism and biosynthesis of protective compounds such as ectoine and were upregulated in fish fed diets supplemented with algae. This study shows the beneficial potential of Gracilaria as a functional ingredient through the modulation of the complex microbial network towards fish health improvement.

The vertebrate gastrointestinal tract (GIT) is a crowded and complex ecosystem inhabited by microbial communities of bacteria, fungi and archaea that over the last decades has demonstrated to be particularly important in the health and welfare of the hosts 1,2 .
Healthy conditions contribute to a balanced and diversified gut microbiota, preventing its dysregulation or dysbiosis, and promoting beneficial symbiotic interaction with the host 3 . In this symbiosis, while the host provides a good environment and nutrient supply, gut microbiota plays a critical role in nutrient digestion and absorption 4,5 , appetite regulation 6 , immune response 7 , protection from pathogenic microorganisms 8 and gene expression regulation 9 . However, gut microbiomes are shaped by several factors, including trophic level 10 , environmental conditions 11 and feed source 12,13 .
Bacteria communities´ density and composition vary along the GIT depending on the physical and chemical conditions 14,15 . Generally, bacterial density increases progressively along with the GIT, being the intestine the region with high alpha-diversity indices 16,17 . In addition, some bacteria species are present in all GIT, varying their abundance, while others are specific to some GIT regions 18,19 . These variations are related to pH 19 , protease activity 20 , amount and type of available nutrients 18 , and adhesion capacity of bacteria groups to the epithelial cells or mucus 21 among other factors. It is important to understand the dynamics of the microbial community in the gut in order to relate with the metabolic and physiological impacts of nutritional and health alterations in fish 12,22 .
In livestock productions such as aquaculture, a balanced microbial community gains particular importance 23 due to the captive breeding and rearing conditions 24 . High stocking densities and high levels of stress can lead to disease spreading/outbreaks 25 . European seabass (Dicentrarchus labrax) is one of the most relevant produced species in southern Europe, and its production is strongly affected by bacterial diseases, mainly photobacteriosis, vibriosis and tenacibaculosis 26 . Regarding the photobacteriosis, the responsible pathogenic agent is Photobacterium damselae subsp. piscicida. This bacterium is considered by worldwide fish farmers as one of the most dangerous microorganisms due to its ubiquitous distribution, high mortality rate and large fish species spectrum 27

Results
Fish performance and mortality. Seabass growth performance, feed conversion ratio and overall health improvement have been previously reported 49 . Briefly, weight gained and feed conversion ratio (FCR) by the end of the trial were not significantly different among groups, however fish fed diet with 5% seaweed inclusion tend to gain more weight and to convert better. There was no mortality throughout the trial all fish presented normal behaviour and reaction to feeding moments.
16S rRNA amplicon sequencing output and microbial community analysis. A total of 5,177,608 raw reads were quality filtered and merged into 1,742,241 sequences (Supplementary Table 1) in an average of 35.4% non-chimeric high-quality sequences. From this, a total of 3315 ASVs were clustered, however, only 571 had more than 2 counts and were used in further downstream analysis. After rarefaction all samples reached a plateau (Fig. 1A) of observed features, indicating proper sequencing depth. Microbial communities' richness (i.e., Chao1 index) was not modulated by diets ( Fig. 1B; P > 0.05) and neither was the Shannon diversity index ( Fig. 1C; P > 0.05). On the other hand, the gut community's beta diversity based on the Bray-Curtis dissimilarities indicated that in the anterior intestine the structure of the gut microbial communities is modulated by the inclusion of Gracilaria gracilis biomass or extract ( Fig. 2A; P < 0.02). However, the bacterial communities in the posterior intestine do not seem to be modulated by diets (Fig. 2B).
The most dominant phylum in both intestine sections was Proteobacteria, followed by Actinobacteria, Firmicutes and OD1 (Fig. 3). In the anterior intestine, a clear reduction in Proteobacteria abundance is observed when fish were fed algae or extract supplemented diet, whereas Actinobacteria and OD1 increase, except in the group fed the extract for the latter. In the posterior intestine, changes are not so evident, but abundance modulation was also observed. Regarding the order, abundances were modulated differently between intestine sections. In the anterior intestine Actinomycetales and Sphingomonadales abundance increased in group ALGAE2.5 and EXTRACT, whereas Bacillales abundance increased only in the EXTRACT group. However, Vibrionales abundance was strongly reduced in all groups with ALGAE or Algae EXTRACT dietary inclusion. On the other hand, in the posterior section, the major differences were observed in the abundance of Enterobacteriales which was reduced in all Algae related groups. In this intestinal section, Vibrionales abundance was higher in the group fed with a diet supplemented with the EXTRACT.
The most abundant and differentially modulated genus were Sphingomonas, followed by Photobacterium, Staphylococcus and Vibrio (Fig. 4). In the first three, abundances in the control group (CTRL) gut community Figure 1. Rarefaction curve (A), chao1 richness index (B) and Shannon diversity index (C) of European seabass gut microbial communities for anterior and posterior intestine. Rarefaction curve indicates number of observed features for anterior (blue) and posterior (orange) intestine depending on sequencing depth. Fish were fed a basal diet with no supplement (CTRL) or supplemented with Gracilaria gracilis powdered biomass at 2.5% (ALGAE2.5), at 5% (ALGAE5) or with the seaweed extract at 0.35% inclusion rate (EXTRACT). No significance differences were observed; dots (in A) and bars (in B and C) indicate mean value of the group and error lines indicate ± SD; AI and PI stand for anterior and posterior intestine respectively (Kruskall-Wallis, P > 0.05). www.nature.com/scientificreports/  www.nature.com/scientificreports/  www.nature.com/scientificreports/ are higher, especially in the case of Photobacterium. However, the highest abundance of Vibrio was observed in the communities of fish fed diet supplemented with seaweed EXTRACT. Correlation network highlighted several potential relationships between genera based on their abundance in each group (Fig. 5). In anterior intestine (Fig. 5A) it is noticeable a central cluster of genera positively correlated and with higher abundance in ALGAE2.5 group, and this includes the genera Enterococcus, Brachybacterium, Stenotrophomonas, Bradyrhizobium, Brevibacterium and Dietzia. Photobacterium genus highlights with higher abundance in CTRL group and negatively correlated with Propionibacterium, and a smaller cluster is visible with genera with higher abundance in the EXTRACT group that includes the Enhydrobacter, Cloacibacterium, Paracoccus, Coxiella and Dermacoccus, all positively correlated among them. In the posterior section, network is more diffused, and a minor cluster is evidenced where genera with higher abundance in CTRL group are positively correlated, and this includes Cloacibacterium, Streptococcus, Haemophilus, Dermacoccus and Kokuria, whereas another smaller cluster evidence genera with higher abundance in EXTRACT group such as Providencia, Finegoldia, Brevibacterium, Coxiella and Gluconacetoba.

Scientific Reports
Random Forest analysis unravelled the most important genera within the groups' intestinal microbial communities. This assessment is based not on the abundance value but on the magnitude of the modulation of the abundance depending on treatments. Thus, in both intestinal sections the genus Photobacterium was highlighted as the most important feature for these communities (Fig. 6A,B). This genus abundance was significantly higher in CTRL group with a negative modulation exerted by diets supplemented with ALGAE biomass and EXTRACT (P < 0.03 and P < 0.0001 in anterior and posterior intestine respectively). In the anterior intestine community, Staphylococcus and Sphingomonas also presented a relevant role (Fig. 6A) mainly due to their abundance increase in the group fed with EXTRACT, whereas in posterior intestine community Dermacoccus, Anaerobacillus and Staphylococcus were the most relevant features of the community (Fig. 6B). Here only Dermacoccus had a reduction of abundance when fish were fed algae related diets, since the latter two genera presented higher abundances when fish were fed 5% seaweed supplemented diets.
Microbiome functional prediction. Microbiome functional profile prediction was performed based on metabolic pathways, and a hierarchical clustering revealed that the microbiomes functions of fish fed CTRL or ALGAE related diets tend to be different, mainly in the anterior intestine (Supplementary Fig. 2A). Here, communities of fish fed with CTRL diet cluster closely (orange samples in heatmap), and except for two samples of fish fed algae supplemented diet all algae related samples clustered (blue and green samples on heatmap) together (regardless if supplementation was with seaweed biomass or extract). However, in the posterior intes- Figure 4. Relative abundance of most abundant genera in European seabass gut microbial communities when fed the experimental diets, and these were a basal diet with no supplement (CTRL) or supplemented with Gracilaria gracilis powdered biomass at 2.5% (ALGAE2.5), at 5% (ALGAE5) or with the seaweed extract at 0.35% inclusion rate (EXTRACT). Bars and errors indicate mean ± SD of differentially abundant genera (ANOVA, P < 0.05). www.nature.com/scientificreports/ www.nature.com/scientificreports/  www.nature.com/scientificreports/ tine ( Supplementary Fig. 2B) this differentiation is not as evident, and samples present more similar functional patterns (samples with mixed clustering). When evaluating the metabolic pathways abundance (i.e., predicted enrichment) in the anterior intestine ( Fig. 7A) 16 out of 17 pathways are enhanced in the microbiome of fish fed either 2.5% or 5% ALGAE biomass supplemented diets. The only pathway that is less enriched in the ALGAE group is the super pathway of l-alanine biosynthesis, whereas several present strong enhancement in this group (e.g. heterolactic fermentation, Bifidobacterium shunt, cob (II)yrinate a,c-diamide biosynthesis II, ectoine biosynthesis, among others). In the posterior intestine (Fig. 7B) less predictive enrichment was observed, however, 12 out of 13 differentially enriched pathways were enhanced in ALGAE groups compared to CTRL. Here we highlight l-isoleucine biosynthesis V, succinate fermentation to butanoate (not the highest mean proportion but with high significance), the pyruvate fermentation to acetate and lactate II and the 3-phenylpropanoate degradation.

Discussion
Seaweed's potential for health management has been highlighted over the last years. The antioxidant, immunostimulant, and overall health-enhancing effects exerted by seaweeds have been pointed to as the main reasons for its health-promoting qualities 61,62 . Although a large amount of the available reports refer to the effects on www.nature.com/scientificreports/ human health, several studies have shown their potential also in fish 49,59,[63][64][65] . Members of Gracilaria genus have been investigated as infeed ingredients to improve health in fish, and despite some differences, overall, the inclusion of this seaweed biomass has a positive effect on fish health. By using 16S rRNA amplicon sequencing we were able to identify more than 150 taxonomical features, and rarefaction indicated that sequencing was deep enough to identify all features in samples by reaching a plateau for all samples. Interestingly, although the total number of observed features was not different between the anterior and posterior intestine, we identified a different composition of the microbial communities from both compartments ( Supplementary Fig. 1) based on the Bray-Curtis dissimilarities. Since differences in microbiome composition and functionality between fish intestinal compartments have been identified in other studies 60,66 , analyses were conducted separately for both sections. We observed that neither seaweed biomass nor its extract inclusion in diet changed community richness or diversity. These results differ from the reported by Rico et al. 67 , where higher richness of seabream intestinal communities was observed after being fed with diets including 15% Gracilaria or Ulva biomass. These differences can be attributed to the practiced inclusion rates. However, the modulatory effects of algae dietary inclusion in intestine microbial communities' diversity have been in some cases contradictory. When considering macroalgae, no differences in diversity in seabream intestines were found by Abdala-Díaz et al. 68 after 30 days of feeding with Ulva rigida at 25% inclusion. On the other hand, Tapia-Paniagua et al. 60 found higher diversity in Senegalese sole anterior intestine microbiota after adding 5% Ulva ohnoi to the diet for 45 days. However, compared with other dietary supplements with potential as fish health modulators (e.g., probiotics, prebiotics), the algae effect on the gut has received limited attention with only a few published studies. Nevertheless, reported results allow us to infer that microbiome modulation by algae is highly dependent on host and algae species, inclusion rate, and duration of feeding among other factors, reinforcing the complexity of the gut microbial community and its relationships.
In terms of composition, Proteobacteria, Firmicutes and Actinobacteria have been reported as the most common phyla of intestinal bacteria in marine fish. The microbial community of seabass evaluated in this study was in line with the previously described in both intestine sections in any diet enrichment 13,23,67,[69][70][71] . At the phylum level, the results showed that the inclusion of Gracillaria gracilis in the fish diet induced an evident change in the anterior intestine, decreasing Proteobacteria abundance, and this is in line with other studies that used Ulva ohnoi 60 . On the other hand, Actinobacteria abundance increased in the anterior intestine of fish fed with algae or algae extract. Other studies have shown different modulations with an unchanged or decreased abundance of this phylum and interestingly, the use of Ulva ohnoi or Gracilaria sp. increased Tenericutes and Firmicutes abundance, respectively 60,71 . The strong reduction of Vibrionales abundance observed in all groups with algae or algae extract is in line with other studies that showed less abundance of Vibrio and Photobacterium genera, both members of Vibrionales order, when fish were fed with algae supplemented diets 23,60,67 . Members of these genera are Gram-negative bacteria that gained notoriety partly due to the pathogenesis of some of its species. In particular, microorganisms from Photobacterium genus can be isolated from various marine surfaces and environments, including other organisms, with which they establish interactions that may be negative to the host 72 . Photobacterium damselae, is a fish pathogen with two subspecies, P. damselae subsp. damselae and P. damselae subsp. piscicida (Phdp) which is the infectious agent of pasteurellosis in fish. Due to its low specificity and high mortality rates, is liable for huge economic impacts in the industry on a global scale 28 . In the present study, the inclusion of Gracilaria as a functional feed additive influenced gut microbial composition, and one of the observed modulations was on the abundance of genus Photobacterium. In fish fed diets with Gracilaria inclusion, but not its extract, the relative abundance of Photobacterium was reduced, limiting its role as a permanent and latent member of the intestinal microbiota. This modulation was highlighted as the most relevant taxonomical modulation in the gut community by the random forest analysis, unravelling a probable relationship between in-feed administration of Gracilaria and Photobacterium species abundance. Although this requires further validation, this result is in line with the ones described by O'Sullivan et al. 73 , which suggested that some macroalgae polysaccharides (i.e., agar and carrageenan) may have inhibitory effects towards some microorganisms. When feeding European seabass with diets supplemented with 5% Gracilaria gracilis extract, Peixoto et al. 65 found an increase in fish resistance to infection with Photobacterium damselae subsp. piscicida, which was explained by a higher antioxidant and immune response. Interestingly, Passos et al. 59 also found an improved resistance against the same pathogen in gilthead seabream fed with diets similar to the ones used in this study. Since the overall physiological analysis did not indicate a direct cause of the resistance (i.e., immunostimulation), it was suggested gut microbiota could have played a role in this feature. A possible antagonistic relationship between the seaweed and this pathogenic bacterium was evaluated by Passos et al. 49 with non-significant results. However, in this study we identified several clusters of correlations between some genera that indicate the bacteria relationships that are modulated by the algae-supplemented diet in both intestinal sections, but not by the algae extract.
In the anterior intestine, Photobacterium presented a possible antagonistic (i.e., negative correlation) relation with Propionibacterium, a microorganism that produces propionic acid, a short-chain fatty acid with antimicrobial properties 74 . There was an evident cluster where the genera Burkholderia and Streptococcus were correlated due to their abundance reduction in fish fed the diets with algae inclusion. Both genera include pathogenic species 75,76 but for example, Burkholderia has also been linked to higher health performance in salmon 15,75 whereas Streptococcus have been highlighted as probiotics in shrimp 77 . Another cluster of positively correlated genera was evidenced due to their higher abundance in fish fed the diet with 2.5% algae inclusion. The correlated genera, Brachybacterium, Enterococcus, Stenotrophomonas, Brevibacterium, Dietzia and Bradyrhizobium are naturally occurring taxa with different properties among them such as lipase production 78 , regulators of nitrogen and sulfur cycle 79 or even opportunistic behavior 80,81 , and might be playing a relevant role in the modulation of the seabass intestinal microbiota community co-occurring in a well-orchestrated manner. Changes of these bacteria abundances and associations might have consequences for the microbiome functioning and its relationship with the host, in this case, promoting better health status and disease resistance 49 www.nature.com/scientificreports/ assessment of fish from this trial, Passos et al. 49 highlighted some histological alterations in the gut. Interestingly, in the groups where microbiota modulation was more noticed (ALGAE2.5 and ALGAE5) histological changes did not show a coherent pattern. When fish were fed ALGAE5 diet villus length were lower, whereas when fed ALGAE2.5 villus width increased. However, all algae or extract fed group had an increase in the number of goblet cells indicating an effect in the mucosa. Although in our study we found no significant correlation between microbial abundance and goblet cell quantification, this was probably due to limited sample number (i.e., replicates average was used since histology and microbial analysis were not obtained from same fish), and it should be targeted in a future study to understand the relation with microbial community. Indeed, to deepen our knowledge on the extent of the observed microbial modulation, a microbiome functional prediction was performed and highlighted a possible modulation of several metabolic pathways in the microbiome of fish fed algae supplemented diet. For instance, an increase in heterolactic fermentation was predicted when including Gracilaria biomass in the diets, and this is probably linked with lactic acid bacteria (such as from the genus Enterococcus) activity and is likely to have beneficial outputs for the host intestinal health 82 . An increase in Bifidobacterium Shunt pathway was also predicted and it might be responsible for an increase of acetate and lactate in the gut lumen. These compounds acidify the intestinal lumen preventing the growth of harmful bacteria, and also serve as an energy source for intestinal epithelial cells 83 . Although the levels of these compounds were not assessed in this study, a further validation considering these assessments would confirm this possible beneficial effect G. gracilis dietary supplementation. Another interesting example is the ectoine biosynthesis pathway that was predicted to increase in the microbiome of fish fed algae supplemented diets. Ectoine is a natural compound found in higher concentrations in halophilic microorganisms and acts as a compatible solute for the survival of osmotic stress 84 . Its commercial form is issued in nutraceuticals as an enzyme stabilizer and cell protector for skin, and a similar role in the gut should be further investigated. It is worth mentioning that in the posterior intestine the modulation observed on the Photobacterium genus abundance was even more noticeable and here it was found with a positive correlation with members of the genera Rubritalea and Gardnerella, whereas the genera Anaerobacillus, Micrococcus and Rubribacter seem to have an antagonistic relationship with the pathogenic genus members. These interactions are yet to be studied and require further confirmation as well as the consequences for the microbiome functioning. However, it is worth mentioning that in the posterior intestine, it was predicted an in the pyruvate fermentation to acetate and lactate II and succinate fermentation to butanoate pathways in the microbiome of fish fed algae supplemented groups. Although these results are yet to be validated with target analysis, an increase in these pathways would lead to an increase in the production of short-chain fatty acids (e.g., butanoate) as well as acetate and lactate. If confirmed this will have beneficial consequences for the gut epithelium. More, understanding how dietary supplements modulate gut microbial interactions and their cross-talk with the host under different functional nutrition scenarios is of utmost importance nowadays. Correlating that information with fish physiological output will allow to move forward in fish intestinal health management in aquaculture.
In conclusion, supplementing seabass diets with Gracilaria gracilis biomass at 2.5% and 5% has an impact on gut microbiome composition. The diet did not alter the diversity and richness of the communities, however, alteration of the abundance patterns of some taxonomical groups was observed and at the genus level, several taxa presented correlation patterns that suggest a possible mutualistic/antagonistic coexistence. The modulation observed in the abundance of members of Photobacterium genus was the most relevant alteration exerted by the diets, with abundance reduction of the pathogen to undetectable levels.

Methods
Ethics statement. The current study was conducted complying with the ARRIVE guidelines and according to the European Directive 2010/63/EU. All experimental procedures were approved by DGAV (Portuguese Veterinary Authority) under the license 0421/000/000/2019 and by the Animal Welfare Committee of the Polytechnic Institute of Leiria.
Algae collection, processing, and experimental diets. Algal biomass of Gracilaria gracilis, harvested from the Portuguese west coast was brought to Cetemares facility (MARE-Polytechnic of Leiria, Peniche, Portugal). All contaminants were removed and G. gracilis was thoroughly washed with seawater. Procedures were then followed to obtain dry algae powder and algal extract from the clean seaweed biomass. Algae powder was produced by drying the seaweed at 25 °C until constant weight and then grinding it to dust. The samples were stored at − 20 °C until use. The algal extract was prepared by drying the algal biomass at 25 °C until constant weight, grinding it into particles smaller than 200 µm, extracting twice with distilled water in a proportion of 1:10 (m:v) and then extracting in absolute ethanol in a proportion of 1:10 (m:v), all extractions were performed in a magnetic stirrer at room temperature for 30 min. The crude extract obtained was then filtered through a paper filter (Whatman nº4) and evaporated in a rotary evaporator at 40 °C. The extracts were frozen at − 20 °C until further processing. A specialized company (Sparos Olhão, Portugal) mixed the algal extract and powder in standard aquafeed 49 , considering an algal extract concentration of 0.35% (EXTRACT) and two algal powder concentrations, 2.5% (ALGAE2.5) and 5% (ALGAE5). Dry ingredients were mixed in a double-helix mixer (model RM90, Mainca, Barcelona, Spain) and ground (below 200 µm) in a micropulverizer hammer mill (model SH1, Hosokawa-Alpine, Augsburg, Germany). Subsequently, the oils were added to the mixtures, which were humidified with water and agglomerated by a low-shear and low-temperature extrusion process (Italplast West Heidelberg, VIC, Australia). Extruded pellets were dried in a vibrating fluid bed dryer (model DR100, TGC Extrusion, Roullet-Saint-Estèphe, France). Diets were packed in sealed plastic buckets and shipped to the experimental facilities (MARE-Polytechnic of Leiria, Peniche, Portugal) where they were stored at room temperature in a cool and aerated emplacement. Proximal composition of the diets was not different and was previously www.nature.com/scientificreports/ reported by Passos et al. 49 and ash content ranged between 7.4 and 8.1%, and in relation to dry matter protein ranged between 52.2 and 53.5%, fat between 13.5 and 13.9% and energy between 21.7 and 21.9 (KJ g −1 )/%DM.
Fish and rearing conditions. The trial was performed at the Aquaculture Laboratory of MARE-Polytechnic Institute of Leiria (Peniche, Portugal) and al procedures were previously approved by the Ethical committee and were previously reported by Passos et al. 49 for the first part of the study. European seabass (Dicentrarchus labrax) juveniles (17.49 ± 6.07 g; mean ± SD) were obtained from Estação Piloto de Piscicultura de Olhão (Instituto Português do Mar e da Atmosfera, I.P.) and were acclimated to the laboratory facilities for two weeks. After the quarantine period, fish were randomly distributed into 12 aquaria (60 L, 5.83 ± 0.31 kg m −3 ), connected to four recirculating aquaculture systems (RAS). The feeding trial lasted for 47 days, and water parameters such as temperature (19.93 ± 0.54 °C), salinity (31.53 ± 0.50), pH (8.48 ± 0.31) and O 2 (87.87 ± 3.36%) (mean ± SD) were monitored daily. Water ammonium and nitrite levels in the tanks were kept below limits (< 0.05 mg L −1 and < 0.5 mg L −1 , respectively). The fish from each treatment (4 diets and standard feed as control) were handfed, to apparent satiation, twice a day (9 a.m. and 4 p.m.). No mortality occurred during the experimental trial.
Sampling procedure. By the end of the feeding period, 3 fish per replicate were sampled. The fish were anaesthetized with 2-phenoxyethanol (0.5 mL L −1 ) and euthanized by anaesthetic overdose and confirmation by severing the vertebral spine in the immediate post cranial region. Fish abdominal skin was washed with 70% ethanol and all the sampling procedures thereafter were performed under aseptic conditions. As the fish were put through fasting on the last day of the feeding period, the intestinal tracts were empty. Anterior and posterior intestine sections of 2 cm were collected separately, frozen in liquid nitrogen and stored at − 80 °C until processing.
DNA extraction and 16S rRNA sequencing. Total genomic DNA was extracted from the anterior intestine and posterior intestine samples using a sterile scalpel to scrape internal contents and mucosa, applying the QIAamp Fast DNA Stool Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions 12 . The integrity of isolated DNA was evaluated by agarose gel electrophoresis and DNA was quantified in a Qubit fluorometer (Thermofisher Scientific, Waltham, USA). Before library preparation, a PCR was performed to ensure the presence of bacterial DNA using universal primers for 16S rRNA (341F-5′-CCT ACG GGNGGC WGC AG-3′ and 785R: 5′-GAC TAC HVGGG TAT CTA ATC C-3′) 85  Bioinformatic analysis and statistics. Sequences were demultiplexed by the sequencing provider inhouse software, and microbiome bioinformatics analysis was performed with QIIME 2 2020.8 86 . Pairedend raw sequences were filtered for quality, merged and chimeras removed by denoising with DADA2 87 . The obtained amplicon sequence variants (ASVs) were aligned with mafft 88 and phylogeny was constructed with fasttree2 89 . Features appearing in only one sample were considered artefacts and were removed. Taxonomy was assigned to ASVs using the q2-feature-classifier 90 using as database the Greengenes 13_8 99% OTUS reference sequences 91 . All sequences assigned for chloroplast or mitochondria were excluded and all results were separated by tissue (i.e., anterior and posterior intestine). For diversity analysis, samples were rarefied (i.e., randomly subsampled to the smallest library without replacement) to 29,959 sequences per sample. Here, alpha-diversity metrics (observed features, Chao1 and Shannon index) and beta-diversity metrics (Bray-Curtis dissimilarity), and Principal Coordinate Analysis (PCoA) were estimated using the q2-diversity script. Differences between the group's alpha-diversity metrics were assessed using QIIME2 significance tests resulting in an evaluation with the Kruskal-Wallis test, whereas the group's Bray-Curtis dissimilarities significance were tested by PERMANOVA which uses the dissimilarities between samples of the same group and compares them to the distances between groups 92,93 . Differential abundance analysis was performed on non-rarefied data; however, a cumulative sum scaling (CSS) was performed to normalize data. This method accounts for heteroskedasticity of feature variance across values and controls the false discovery ratio in data 94,95 therefore it is preferable to classical total sum scaling. A correlation network analysis was performed to identify possible interactions between microorganisms. Highlighting these interactions can provide valuable inputs on the microbiome function and if a specific diet promotes different interactions between taxa. SparCC 96 correlation method was applied, considering 100 permutations, and retaining features with a correlation coefficient higher than 0.4 and respective P-value < 0.05. SparCC uses a log-ratio transformation and identifies taxa pairs different from background correlations by performing multiple iterations. Correlation analysis was performed based on FastSpar implementation available from the MicrobiomeAnalyst 97 . To identify microbial taxa that differentiate between groups the Random Forest algorithm was applied. This is a supervised machine-learning algorithm that identifies non-linear relationships by constructing multiple decision trees using a randomly selected subset of the data, allowing classification and selection of important features 98 . The random forest model was created using 500 trees. www.nature.com/scientificreports/ Microbiome functional profiling prediction was performed with the Phylogenetic Investigation of Communities by Reconstruction of Unobserved States, named PICRUSt2 99 . This method is based on the idea that phylogenetically related organisms are more likely to have similar gene contents, and the algorithm uses several gene family databases. In this study, the functional profiles of the bacterial communities were predicted using the PICRUSt2 from the MetaCyc pathways database. The differences in functional profiling for the microbial communities in the anterior and posterior intestines of European seabass fed with different diets were characterized using Statistical Analysis of Metagenomics Profiles 100 . The significance level was always set to 0.05 considering a corrected P-value for false discovery ratio (FDR-corrected).

Data availability
The datasets generated and analyzed during the current study are available in the NCBI-SRA repository under the BioProject PRJNA781597. www.nature.com/scientificreports/